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Exact gauge invariant mass dependence of a s through two loops 
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A physically defined QCD coupling parameter naturally incorporates massive quark flavor thresholds in a gauge 
invariant, renormalization scale independent and analytical way. In this paper we summarize recent results for 
the finite-mass fermionic corrections to the heavy quark potential through two loops leading to the numerical 
solution of the physical and mass dependent Gell-Mann Low function. The decoupling-, massless- and Abelian- 
limits are reproduced and an analytical fitting function is obtained in the V-scheme. Thus the gauge invariant 
mass dependence of uv is now known through two loops. Possible applications in lattice analyses, heavy quark 
physics and effective charges are briefly discussed. 
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1. Introduction 

Quark flavor thresholds in QCD are commonly 
treated within effective descriptions in MS-like 
coupling definitions by imposing matching condi- 
tions at the quark thresholds Jj]||. Thus quarks 
are considered infinitely heavy below and mass- 
less above m q and the coupling is non-analytic 
at the thresholds. Real mass effects need to be 
calculated separately as higher twist effects in 
the small and large mass limits. For the in- 
termediate range an all orders resummation of 
these expansions is necessary. In this paper 
we summarize recent results presented in Refs. 
|U^] based on a physical coupling definition ob- 
tained from the static quark-antiquark potential 



which natu- 



@, U(Q 2 ,m 2 ) = -47rCj 

rally incorporates massive quarks and where the 



ay(Q 2 ,m 2 ) 
Q 2 



scale Q 



-q = q is identified with exchanged 



momentum between the heavy sources. A tech- 
nical complication is that the massive Gell-Mann 
Low function can only be solved numerically due 
to the complexity of the obtained results and 
that it is scheme dependent already at one loop. 
The latter point can be ameliorated by express- 
ing other physical charges through ay and using 
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the conformal ansatz Q. We begin in the next 
section by reviewing the two-loop corrections in- 
cluding massive quarks to V(Q 2 ,m 2 ) and then 
discuss the solutions to the massive renormaliza- 
tion group equations. Finally we briefly outline 
possible applications. 



2. Two loop corrections 

The results obtained in Ref. Q express the 
physical charge ay in the MS-scheme, which is 
used as a calculational tool, in the following way: 



a v (Q,m) = a MS (/i) + rn(^i), 
+ U2 (Q,m( M ), A i)^4^ 



(1) 



where v-i contains the diagrams of Fig. |l| and the 
MS-counterterms displayed in Fig. g. A strong 
check of the results in Ref. || is given by the suc- 
cessful reproduction of the fermionic gluon wave 
function renormalization constant (RC) and the 
locality of all other RC's as these are mass inde- 
pendent in minimally subtracted schemes. 

For the heavier quark masses m c , and nit 
the pole-mass definition is suitable and allows for 
a straightforward Abelian limit as well as the 
renormalization scale independence of the Gell- 
Mann Low function below. The next-to-leading 
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Figure 1. The massive fermionic corrections to 
the heavy quark potential through two loops in 
the Feynman gauge. The straight ladder diagram 
does not contribute as it is already contained in 
the iteration of lower order amplitudes. The mid- 
dle line contains IR-divergent diagrams, however, 
their sum is IR-finite. Contributions proportional 
to Cf and Ca are separately gauge invariant for 
ay. After inclusion of the counterterms in Fig. ^ 
the correct massless limit given in Refs. ||[7| is 
obtained. Details can be found in Ref. 101. 



Figure 2. The two loop counterterms correspond- 
ing to the diagrams in Fig. [|. Adding these 
contributions to the original graphs removes all 
non-local functions from the occurring pole terms. 
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The only exception are — terms in the two point 
functions which only cancel in the sum of all two 
point diagrams. The fact that the tadpole dia- 
gram has no counterterm is already indicative of 
this cancellation. 



order relation between the MS mass m{\i) and the 
pole mass m is given by S : 



m(/i) 



1 - C F 1 + - log — 

7r \ I m 



- [ 7 - log(47r)] 



(2) 



where 7 is the Euler constant. Inserting Eq. (|2|) 
into Eq. (|l|) gives at next- to-next-to-leading order 

,«ms(a*) . 



av(Q,m) 



1 + vi(Q,m,n)- 



[v 2 (Q, m, fi) + A m (Q, to, fx)] 



(3) 
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other words we will write 
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Figure 3. The numerical results for the gauge- 
invariant Npy in QED (open circles) and QCD 
(triangles) with the best \ 2 fits of Eqs. ( |l0| ) and 
(||) superimposed respectively. The dashed line 
shows the one-loop Np V function. For com- 
parison we also show the gauge dependent two- 
loop result obtained in MOM schemes (dash-dot) 
|lO|JTl| l . At large — the theory becomes effectively 
massless, and both schemes agree as expected. 
The figure also illustrates the decoupling of heavy 
quarks at small — . 



where A m (Q, to, denotes the contribution aris- 
ing from vi when changing from the MS mass to 
the pole mass. 

3. Numerical solutions of the Gell-Mann 
Low function 

The Gell-Mann Low function g| for the V- 
scheme is defined as the total logarithmic deriva- 
tive of the effective charge with respect to the 
physical momentum transfer scale Q: 



day(Q, m) 
dlogQ 



i=0 



( + 1 



(4) 



For the massive case all the mass effects will be 
collected into a mass-dependent function Np. In 



4 0) 



11 

y 



(0) 
F,V 



r 2l _ Hjv w 
4 n F - v 



Q 

m 

Q 
m 



(5) 
(0) 



where the subscript V indicates the scheme de- 
pendence of Np V and Np V . 

Taking the derivative of Eq. (||) with respect 
to log Q and re-expanding the result in av(Q, to) 
gives the following equations for the first two co- 
efficients of ^!y. 



\ TO 



(7) 



dlogQ 

d[v 2 (Q, to, n) + A m (Q, to, fj,)] 
dlogQ 

+2v 1 {Q,m,fi) — - — - — 8 
dlogQ 

The argument Q/m indicates that there is no 
renormalization-scale dependence in Eqs. (Q) and 



(|8j). Rather, ipy^ an d ^v' are functions of the ra- 
tio of the physical momentum transfer Q = \J — q 2 
and the pole mass to only. A numerical solution 
based on the MC-integrator VEGAS and numeri- 
cal differentiation gives stable results summarized 
in Fig. |^. In the case of QCD we obtain the fol- 
lowing approximate form for the effective number 
of flavors for a given quark with mass to B: 



N 



(i) 

F,V 



-0.571 + 0.221 



Q 2 \ Q 



m* 



(9) 



1.326-+T 

TO 



0.221^ 

TO 



and for QED 



N 



(i) 

F,V 



1.069 + 0.0133- 



Q 2 \ Q 



,Q 4 



0.402-^r + 0.0133-^r 

TO TO 



(10) 



(1) 

F,V 



The results of our numerical calculation of N j 
in the F-scheme for QCD and QED are shown in 
Fig. ||. The decoupling of heavy quarks becomes 
manifest at small Q/m, and the massless limit is 
attained for large Q/m. 
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Figure 4. The upper diagram displays a compar- 
ison of the Abelian limit of our results (open cir- 
cles) for Npy based on the calculation in Rcf. H 
which was done in the MS-scheme with the well 
known result in the literature jjj] done in the on- 
shell renormalization scheme (solid line). Also 
shown are the gauge invariant non- Abelian con- 
tribution only (oc Ca) (open triangles) as well as 
the sum of all terms in QCD (solid triangles). The 
correct Abelian behavior is a very strong check on 
the results given in Ref. |). The lower dia gram il- 
lustrates the renormalization scale independence 
of the two-loop effective number of flavors Npy 
as a function of the ratio of the physical momen- 
tum transfer Q over the pole mass m. Numeri- 
cal instabilities are visible for small values of ^ 

rn 

and occur because of limited Monte Carlo statis- 
tics (10 7 evaluations for each of the 50 iterations). 
The two fits obtained, which agree within statis- 
tical errors, are shown as a solid and dashed line 
for fj, = m and /i = 0.031m respectively. 




Figure 5. The upper plot shows the sum of the 
effective number of flavors for one (dashed line) 
and two loops in the V-scheme. We use quark 
pole masses with m c = 1.5GeV, mb = 4.5GeV 
and rat = 173GeV. The two loop Np V starts to 
decrease from the fixed starting point Npy = 
3 due to the novel non-Abclian anti screening 
corrections and then increases more rapidly as 
the one loop Npy. Below ITeV, there is no 
regime for which the quark masses can be ne- 
glected. The lower plot displays the scaled i Sf- 
function, — Wy/(a v /w) m tne analytic V-scheme 
ay(Q,rrii) (solid) compared to the ay(Q,G) 
scheme with discrete theta-function treatment of 
flavor thresholds with continuous matching at 
Q = m (dashed). 
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We can also apply the same fitting procedure to 
the dependence of the one-loop effective Np V « 
— t ■ Fig- demonstrates that the new non- 

1+5-19^ 6 Q 

Abelian contributions (~ Ca) are responsible for 
the negative Np V at intermediate Q/m due to 
anti-screening. The Abelian corrections on the 
other hand are larger than 1 in this regime and 
agree with the literature results obtained in the 
on-shell scheme ||. In addition the lower graph 
of Fig. H demonstrates the renormalization scale 
independence of the solution to the Gell-Mann 
Low function. 

Fig. |^ demonstrates the smoothness and an- 
alyticity of the renormalization group solutions 
and compares the massive results with the mass- 
less ones including one-loop matching at the two 
loop order. The figure demonstrates that there is 
really no regime below 1 TeV where quarks can 
be considered massless for running coupling ef- 
fects in the V-scheme. 

4. Conclusions 

In summary, we have presented the gauge in- 
variant mass dependence of a s through two loops 
in the physically motivated V-scheme. The result 
was shown to posses the correct massless limit 
and gives automatic decoupling of heavy quarks. 
In addition the correct Abelian limit is repro- 
duced and the renormalization scale independent 
results can be parameterized by a simple analyt- 
ical fitting function. Non-Abelian anti-screening 
effects lead to a negative number of flavors for 
intermediate energies at the two loop level. 

Massive renormalization group solutions are 
scheme dependent already at one-loop, however, 
the mass dependence of ay can be transferred 
to other physical charges through commensurate 
scale relations In Ref. Q this was done for 
the non-singlet hadronic width of the Z-boson and 
compared with the MS'-scheme higher twist cor- 
rections. For perturbative energies a persistent 
~1% deviation was observed which characterizes 
the residual scheme dependence of the two loop 
predictions for this observable. In a similar way, 
all other effective charges can be described. 

Other possible applications include the effect 
of a massive charm on the bottom mass deter- 



mination. Here massive charm corrections to the 
potential and in the running coupling could po- 
tentially lead to a shift in the bottom mass of 
O (m c cxj.(mb)) and thus would need to be in- 
cluded into a proper analysis. Also top quark 
physics at the NLC could provide fruitful ground 
for a V-scheme analysis. 

A further interesting comparison could be 
performed with lattice analyses investigating 
the transition region of perturbative and non- 
perturbative regimes. For this purpose the 
Fourier transform ay(l/r) = J ay(q 2 , m 2 ) 



= „2 av(l/r) 



must be obtained. 
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